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Abstract 



We consider the properties of the diffusion controlled reaction ^ + S — > in the 
steady state, where fixed currents of A and B particles are maintained at opposite 
edges of the system. Using renormalisation group methods, we explicitly calculate 
the asymptotic forms of the reaction front and particle densities as expansions 
in {JD~^\x\'^^^)~^, where J are the (equal) applied currents, and D the (equal) 
diffusion constants. For the asymptotic densities of the minority species, we find, 
in addition to the expected exponential decay, fluctuation induced power law tails, 
which, for d < 2, have a universal form 74|a;|~'^, where lv = 5 + 0(e), and e = 2 — d. 
A related expansion is derived for the reaction rate profile R, where we find the 
asymptotic power law R ~ B\x\~'^~''^. For d > 2, we find similar power laws with 
LU = d + 3, but with non-universal coefficients. Logarithmic corrections occur in 
d — 2. These results imply that, in the time dependent case, with segregated initial 
conditions, the moments / |a;|'^_R(a;, t)dx fail to satisfy simple scaling for g > a; + 1. 
Finally, it is shown that the fluctuation induced wandering of the position of the 
reaction front centre may be neglected for large enough systems. 



1 Introduction 



Since the initial work of Galfi and Racz [|T|, there has been considerable interest 
in the kinetics of one and two species annihilation, A + A ^, and A + B —>■ ^ 
[1-17, 23-27]. Most analytic and numerical studies have concentrated on the case 
of either homogeneous initial conditions, or initially entirely segregated reactants. 
Ben-Naim and Redner were the first to study the case of a steady state reaction 
interface, maintained by fixed particle currents imposed at opposite edges of the 
system. Their equations for the particle densities a(x, t) and 6(x, t) were 

^ = DV^a - Xab (1) 
dt ^ ^ 

^ = DV% - Xab, (2) 
dt ^ ^ 

with diffusion constant D, reaction rate constant A, and with the boundary condi- 
tions: 

J = -Ddxa\^=-L = -DdJ)\x=-L = -Dd:ca'\x=L - J = -DdJ)\^=L. (3) 



These equations are asymptotically soluble analytically, giving 



(4) 

where 9{x) is the Heaviside step function. The relations w ~ J^^^^, where w is 
the reaction front width, and c ~ J^^^, where c is the particle concentration in 
the reaction zone, are also derived in However, implicit in their formulation is 
the 'mean field' like assumption, (ab) oc {a){b), which will no longer be adequate 
below the critical dimension, due to fluctuations. Cornell and Droz |T2| have given 
an argument for the upper critical dimension of the system (leading to = 2), 
as well as performing numerical simulations. On the basis of these, and mean 
field analysis, they have proposed scaling forms for a,b and the reaction front R, 



which are postulated to be vahd both above and below the critical dimension in 
the scaling limit w —>■ oo (or J 0): 

R=-S[-] a = —A[-) h=—B[-]. 5 

w \wj D \wj D \wj 

In other words, the profiles are characterised by a single length scale which itself 
is suggested in |[T^ to vary as w ~ J~^/^ in = 1, and w ~ J~^/^ for > 2 in 
the scaling limit. Cardy and Lee [EBl have given RG arguments which support 



this conclusion. However, we defer further discussion, especially with regard to the 
presence of multiscaling, until section 6. 

In this paper, we present the results of the first renormalisation group cal- 
culation for the asymptotic properties of the densities and reaction front in the 
steady state, which systematically takes into account the effect of fluctuations in 
the stochastic particle dynamics. Previously, the RG had been used to study the 
late time behaviour of reactions with homogeneous initial conditions (see |TB| and 



references therein). Our calculational framework will bear considerable similarities 
with |T^. The basic plan is to map the microscopic dynamics, in the form of a 
master equation, onto a quantum field theory. This theory is then renormalised 
(for (i < 2) by the introduction of a renormalised coupling, which is shown to have 
a stable fixed point of order e. We then group the Feynman diagrams into sets 
whose sums give a particular order of the renormalised coupling constant. It will 
be demonstrated that this grouping is given by the number of loops. These dia- 
grams may then be evaluated (asymptotically) and the Callan-Symanzik solution 
used, to obtain perturbative expansions for the densities and reaction front. Note 
that for d> 2 no renormalisation is necessary and the diagrams may be evaluated 
directly. 

We now present our results for the asymptotic forms of the densities and reaction 
front profile. It will be shown that at zero loops we find a stretched exponential 



dependence which we include in the following summary, even though we expect 
its effects to be overwhelmed by leading and sublcading power law terms. In 
addition, for o? < 2, we do not rule out the possibility of logarithms in higher order 
terms summing to give a modification to the leading power law given below (which 
results from the straightforward evaluation of the one loop contributions) . So we 
find asymptotically as |x| — > oo for d < 2: 



■{ay 



(J/L')|x||^^~^^| + A,|x|(7-5<^)/i2e-^.l-l^''+^^^^ + A,\x\-'+'^ + ... (6) 
R = X{ab) = A4k|(^'^-5)/^2g-A.|.i«^+i)/2 ^ A,\x\-'+^' + ... (7) 



where 



A -n3787r7/D^i9/i2(4^[!)!^.- , , _ (2c^+l)(2ci+2) ^ 

for d^2: 

(10) 

R = S4(ln |a;|)-7/12|^|3/4g-B.(In|x|)-V2|,|3/2 ^ ^^^^^-7^^ \x\ + . . . (11) 



where 



^1 = 0-3787^^^ B.^liMJ/D))y^ Bs=^-^ (12) 



B, = 0.3787{47ry/''D{J/Dr/'' B, = 



(13) 



and finally, for d > 2: 

{J/D)\x 



(a) 



{\x)} + ^il^r'^'^"""''''''' + C3\x\-'-' + ... (14) 
R = C4|a;|^/S-^^l^l''' + C,\x\-^-' + ... (15) 



where 



^1 = 0-3787|^^^ C, = -{XJlD-f'- (16) 



(A/L')5/12 



^3 = (AJ/D2)-i2-^-'^7r-('^+^)/2(d - l)r <^4 = 0.3787A 

(17) 

Cs = A(A/L>)-2( J/L»)-i2-i-'^7r-('^+i)/'r (^"y^j (^^ - + 3)(rf + 4). (18) 

The layout of this paper is as follows. In section 2, the system is defined using 
a master equation, which is then mapped to a second quantised representation, 
and then to a field theory. In section 3, we present two related field theories and 
derive the form of their Green functions. The renormalisation of the theory is also 
addressed. The calculations for the densities and reaction front are presented in 
section 4, for d < dc, d — dc, and d > dc- The separate problem of the fluctuations in 
the position of the centre of the reaction front is presented in section 5. A discussion 
of these results and comparisons with the available data from simulations are given 
in section 6, where we also argue for the presence of multiscaling in the system. 



2 The Model 

We consider a model where A and B particles are moving diffusively on a hypercubic 
lattice, with lattice constant There is some probability of mutual annihilation 
whenever an A and a B particle meet on a lattice site. In addition, particles of type 
A are added at a constant rate to lattice sites on the hypersurface x — —L, and 
particles of type B arc similarly added to sites at x = L. In other words, opposing 
currents of A and B particles are maintained at opposite edges of the system. The 
two hypersurfaces x — ±L mark the boundaries of the system beyond which the 
particles are not permitted to move. The model is defined by a master equation 
for P{{n, m}, t), the probability of particle configuration {n, m} occurring at time 
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t. Here {n,m} = (ni, n2, . . . , n^r, mi, m2, . . . , m^v), where rii is the occupation 
number of the A particles, and rrii the occupation number of the B particles, at 
the ith lattice site. The appropriate master equation is 

^P({n,m},t) = 
^ + . . ni - 1, Ue + 1, . . . , {m}, t) - niP{{n, m}, t)} 

+ + l)P({n}, . . .nii - l,me + 1, . . . ,t) - miP{{n,m},t)} 

i,e 

+^Yl{i^i + + ■ ■ni + l,mi + 1 . . . ,t) - nimiP{. ..ni.rrii.. . 

i 

+R{P{. . . , n_i - 1, {m}, t) - P(. . . , n_i, {m}, t)} 

+R{P{{n}, mL-l...,t)- P({n}, . . . , t)}, (19) 

where i is summed over lattice sites, and e is summed over the nearest neighbours 
of i. The first, second, and third lines of the equation describe diffusion of the A 
and B particles respectively (with equal diffusion constants D), whilst the fourth 
line describes their annihilation within the system (with rate constant A). The final 
four terms are due to the addition of A and B particles at the edges of the system 
at a rate R, corresponding to the maintenance of steady particle currents. 

The master equation can be mapped to a second quantised form, following a 



standard procedure developed by Doi |Tq] and Peliti |19|, and as described by Lee 



r^. In brief, in terms of the creation/annihilation operators a,a'^,b,b'' which are 
introduced at each lattice site, the time evolution operator for the system is 

H = -^Y.{('l(^e-a,)+bl{K-k)}-XY,{l-albl)aA-R{alL-Mi-l}. (20) 

i,e i 

This can now be mapped onto a path integral, in which a, a, b, b are replaced by 
continuous c-number fields, with action (up to a constant) 
rt ( , . D 



S = Y^iJ dt\ dihi + bibi - — 



di ^{tte - ai) + bi ^{be - bi) 
A(l - dibi)aibi} - ai{t) - kit)) , (21) 
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where aj(t) and bi{t) are due to the projection state (see [0), provided that 

y = ;^K-a_.) = - 



i? = T = 7^K-«-i) = -(6o-6_i) (22) 



at — L, and 



R = j = ^{bo-bL) = ^{ao-aL) (23) 



at +L. Here the sites —L, L are at the edges of the system, with site o being 
immediately outside L or —L. Taking the continuum hmit of this action, we arrive 
at 

S[a, a, b, b,t]= J (^dxd'^-^y j' dt [d{dt - DV^)a + b{dt - DV'^)b 

-Ao(l - db)ab] - a{t) - 6(t)) . (24) 

subject to the conditions 

- J = -DdJ = -Dd^a (25) 

at +L, and 

J = -Dd^a = -Dd:,b (26) 

at —L. Here y are the coordinates for directions perpendicular to the applied 
currents. These conditions may be made explicit in the action by including a pair 
of delta functions: 

S = j {dxd'^-^y y* dt [d{dt - DV^)a + b{dt - DV^)b (27) 
-Ao(l - db)ab - dJ5{x + L) - bJ5{x - L)} - a{t) - 6(t)) . 

If we make the substitutions d = 1 + d and b = 1 + b, then the action becomes (up 
to a constant) 

S = J dxd'^-^ydt [a{dt - D\/^)a + b{dt - DV^)b + \oaab + \obab 

+Xoabab - aJ6{x + L) - bJ6{x - L)] . (28) 



If we integrate over the a and h fields, and neglect the ahab term, we obtain the 
classical (mean field) equations 

{dt - DV'^)a + Aoa6 - J5{x + L) = (29) 
{dt - DV^)b + Xoab - J6{x -L) = 0. (30) 

On the further conditions that no particle annihilation occurs at the edges of the 
system, and that Va = and Vb = outside, integrating the first equation from 
—L — e to —L + e and the second from L — e to L + e in the limit e — > gives the 
required boundary conditions. 

As the diffusion constant exhibits no singular behaviour in the renormalisation 



of the theory, it is convenient to absorb it into a rescaling of time, as in |[T6|. Defining 
t = Dt, A = \oD~^, and J = JD~^, and introducing the fields = |(a + b) and 
= |(a — 6), we have 

S = J dxd'^-^ydt [2^{dt - V^)^ + 2ip{dt - + 2A0(02 - ^2)+ (31) 
A(02 - ^^){(f)^ - - J^[6{x + L)+ 6{x - L)] - J^[6{x + L) - 5{x - L)] . 

Consequently, the new classical equations of the steady state are 

dli^c + ]^J[6{x + L)-5{x-L)]=Q (32) 
dl'Pc - K€ - V^c) + + L) + 5{x - L)] = 0. (33) 

The appropriate solution for ipf. is just —{J/2)x (for |a;| < L), whilst substituting 
4>c = (-^/2)|a;| + u into the second equation gives asymptotically the Airy equation 
for u, as noted in P]. Asymptotically one finds 
0c =^\x\+ 0.3787 (-^y'' |,|-i/v(2/3)(A..-)V.|.|3/. ^ _ _ _ ^ ^34^ 



where the constant was determined numerically. 

So far the quartic term in the action has been neglected, with the result that the 
simple mean field results have been recovered. However, we can take into account 



the non-classical term by including Gaussian noise in the equations for and ijj, 
leading to equations which are exact. This modification can be derived by replac- 
ing the quartic piece in the action by a noise variable, integrating over the noise 
distribution, and demonstrating that this recovers the original term. Observing 
that 

r rfr/<^e^''^e-^^/[^^(<^'-'^')l ~ e^^'^^'"^'), (35) 

J — oo 

and 

J —oo 

where rj^ and t]^ are complex Gaussian noise variables with an appropriate phase, 
we see that the steady state equations may be written as 

dli^ + + L)- S{x -L)] + r]^ = (37) 

dl(f) - A(02 - ^2) + ^J[6{x + L)+ 6{x - L)] +r]^ = 0. (38) 

Clearly we have lost the simple interpretation of the a and b fields as being the 
local densities of A and B particles, as now each of the above equations includes a 
(generally) complex noise term. Nevertheless, we can still interpret {ijj) and (0) as 
being averaged densities, which also satisfy 

dim + + L)~ 5{x - L)] = (39) 

dl{<P) - - ij') + y[5ix + L) + 5ix - L)] = 0. (40) 

The second equation will be used later on to relate a perturbation expansion for 
(0) to one for (0^ — -0^). 

Finally, we give the natural canonical dimensions for the various quantities 
appearing in the action, noting that the coupling becomes dimensionless at the 
postulated value of the critical dimension |T^ : 

[t] = k-^ [a, b] = k'' [a, b] = [A] = k^"" [J] = k''^\ (41) 
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3 Field Theory Formulations 

In what follows it will be convenient to develop two parallel field theories - one given 
by the action already described in (|3T|), and another to be described below, formed 
by writing = 0^ + 0i and ip = ipc + particular, whilst the second theory is 

more useful for calculations, the cancellation of divergences after renormalisation, 
and the identification of leading terms in an expansion in powers of the coupling 
constant, are easier to see in the first theory. 

3.1 Propagators and Vertices 

The propagators for the first of the two theories described above (which we shall 
call Field Theory I) are (from (|3ll) ) 

G^^ikt} = G^4k,i) = h-''\ (42) 

in {k, i) space. In {k, s) space, where a Laplace transform of time has been per- 
formed, we have 

1/2 

G^f{k, s) = G^^{k, = p _^ ^ - (43) 
The vertices are shown in figure 1, where the propagators are represented by 
solid lines, and ip propagators by dotted lines. 

For Field Theory II, we split the and ip fields into their classical and non- 
classical components, which leads to a modified action 

S = J dxd'^-^ydt {20((9f - + 2A0e)0i + 2ij{dt - V^)Vi (44) 
+2A0(0? - ^pf - 2^,^lj,) + A(0' - ^2) (02 + 20,01 + 02 - - 2^,^, - ^i)} , 

where the classical equations have been used to simplify its form somewhat. We 
can now substitute for the exact value of tpc = —{J/2)x and for the functional form 
of the 0c field (from (H)) 

^c = {JVxy/'f[{xjY^H (45) 
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If we also make the rescaling in the action of 

X = {\jf'\ y = {XJf'^y i= {\jf/H, (46) 
then it is transformed to 

S = j dxd^-^ydii.Xjyi \2^{di - + 2/(x))0i + 2^(9,- - 

+ (02 - ((JVA)^/^(i) + 2/(i)0i + + - , (47) 



where 



h{x) = [f{x)r - \x\ (4^ 



which is essentially just the classical profile of the reaction front. The form of the 
propagators is now 

%0 = ^(AJ)'/'G(x,x',y,y',t), (49) 

where 

[di - + 2/(x)]G(x, x\ y, y', i) = 5{x - x')5{y - y')5{t - i'); (50) 

and 

G^,^(^,5) = ^(Ajy/^^, (51) 

in {k, s) space, or 

-{fci+S)l/2|i-a.'| 

G^,j;(x - x', k^, ~s) = {XJf'^ ^- — , (52) 

^i^v , ±, ; V ; 4(A;i + s)V2 ^ ' 

in (x, A;^, s) space, where the perpendicular directions are defined to be those per- 
pendicular to the apphed currents. Unfortunately the equation for G (|50|) is too 
hard to solve exactly, as we do not have an analytic form for 0c- Consequently we 
must rely on the approximation f{x) ~ valid at large in order to make 
the equation tractable. If we also Laplace transform time, and Fourier transform 
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to momentum space for spatial dimensions perpendicular to the applied currents, 
then we obtain 

{s+~kl-dl + \x\) G{x, x', ~k^, s) = 6{x - x'). (53) 

This has the solution (for x' > 0, and accurate for large \x\ i.e. when |x| ^ 
(AJ)-V3): 

' aAi[—x + kj_ + s] when x < 

G{x,x\k±J) = < (3Ai[x + kl + s] +-fBi[x + kl + s] when < x < x' (54) 

, + kj_ + s] when x > x' . 

Considering the boundary conditions at x' (continuity in G and a discontinuity in 

its derivative), we have 

(5Ai{x' + A;^ + 5) + '^Bi{x' + A;^ + 5) = bA%{x' + fc^ + s) (55) 
^A^{x' + A;i + 5) + -iBi!{x' + A;i + 5) - 6Ai!{x' + fc^ + s) = 1. (56) 

These equations can be solved for 7 with the result that 7 = TTAi{x' + k'^ + s). The 
final boundary condition (G — > as x — —00) will (in principle) give a further 
relation between P and 7, as well as specifying a. But to use this condition we need 
to know the behaviour of G in regions near 0, where our asymptotic approximation 
breaks down. Consequently we must rely on numerical solutions, which reveal that 
for our purposes we may neglect the jSAi term in (|55|). Solving for 6, we obtain: 



nAi(x' + k^: + s)Bi(x + k"^, + s) for < x < x' 
G{x,x',~kl,s) = { ^ ^ (57) 

■KBi{x' + kl + s)Ai{x + k'i + s) for x > x' > 0. 



We can now use the asymptotic form of the Airy function [20| to simplify these 
expressions further: 

Hence, for <^ x < x', 

G = i(£' + kl + s)-'/\x + kl+ 5~)-l/4g-(2/3)[(S'+fci+«-)3/2_(,+fc2 +,)3/2]^ ^^g^ 
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with a similar expression for x > x' ^ 0. For x ~ x' ^ 0, we may expand the 
terms inside the exponential, to obtain 

G,^^{x,x'M,~s) = i(Ajy/='(£ + ^i + 5)-^/V(*+^i+^-)^^^l^-^'l (60) 

The important point to notice here is that for x sufficiently close to x', the Green 
function decays only as a power law. 

Finally, we note that each occurrence of a propagator is associated with a factor 
of (XJ)'^^^. If we also extract a factor of (A J)"*-'^"'"^''/^ from each vertex, then we can 
use the vertices shown in figure 2, provided we multiply any given diagram by a 
factor of 

(xj)ip'i-m^^ (61) 

where p is the number of propagators and v is the number of vertices. Again, in 
figure 2, cj) propagators are solid lines and ip propagators are dotted lines. Note the 
simple form of the vertices (h) to (m). 

3.2 Renormalisation 

The renormalisation of the theory proceeds in a similar vein to that described in 
[0 - our field theory differs only in the nature of the boundary conditions. Again 
the only renormalisation required is coupling constant renormalisation, as the set of 
vertices for Field Theory 1 allows no diagrams which dress the propagator. Hence 
we have no field renormalisation and the bare propagators are the full propagators 
for the theory. 

3.2.1 Renormalisation of the Coupling 

The temporally extended vertex function for A + 5 — > is given by the sum of 
diagrams shown in figure 3. This sum may be calculated exactly, as done in |jl7| 
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(remembering extra factors of two resulting from the presence of two different types 
of propagator): 

= l + |Ai?2r(e/2)(. + iF)-/2' ^62) 
where B2 = 2/(87r)'^/^, and e = 2 — d. However, as we are now in the time 
independent state, we take s = 0, leading to 

^^^^ = l + iW(e/2)2^/2fc-^- ^''^ 
The vertex function can now be used to define the renormalised coupling, with 



k = K as the normalisation point (differing from [T^). So we have 



gR = K '\{k)\k=^, go = K'X, (64) 

for the dimensionless renormalised and bare couplings respectively. The f3 function 
is defined by 

f^ign) = f^^gn = -egR + \eglB2T{t/2)T'\ (65) 
and we have a fixed point f3{g'^) = when 

g*j, = {2-'/'B,T{e/2)}-\ (66) 

The fixed point is of order e. Finally, the expansion of go in powers gR remains, as 



m 
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^7o = to + 4 + ■ ■ ■ (67) 
9r 

3.2.2 Callan-Symanzik Equation 

We now write down the renormalisation group equation for {4>i)r (the renormalised 
value of (0i)), expressing its lack of dependence on the normalisation scale: 

f«:^ + /5(to)^') {(pi)R{x,gR,Kj)=0. (68) 
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In addition, dimensional analysis implies 
( d d - d\ 

( K— - X— + {d+l)J^\ qr, k, J) = d{(f)i)R{x, gR, k, J). (69) 

Eliminating the terms involving k, we have 

Xg^-id+ l)JgJ + (^i9R)g^ + d] {<Pl) r{x , Q R, K, J) = 0. (70) 

This can be solved by the method of characteristics, with solution 

(0i)K(a;,fi'ij, J) = (ra)~'^(0i)ij(K"\^ij(K"^),/s;, J(k"^)), (71) 
and associated characteristics 

= /^fe) 9r{x) = gR (72) 

dJ ~_ ~_ 

x— = -(d + l)J J(x) = J. (73) 
ox 

These equations have the exact solutions: 

Jix')^{-) J (74) 

U^')=9Ri^ + ^K^] , (75) 



where in the large |x| limit qr ^ g*^. 

We can make use of the mechanics developed above by first calculating an 
expansion in powers of g^, which can be converted into an expansion in powers of 
gR via (|67|) . Provided that the expansion is non-singular in e, we can relate the gR 
expansion to an e expansion using (0), where for large |x| we can take cjR g*^. 

3.2.3 Tree Diagrams 

At this point we need to identify the leading terms in an expansion in powers of 



g^ - something which can be done in a very similar fashion to 0], using Field 
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Theory I. For the calculation of (0), tree diagrams are of order g^J^'^'^, for integer 
i. Diagrams with j loops will be of order g^J^^'^'K As the addition of loops makes 
the power of g^ higher relative to the power of J, we see that the number of loops 
will give an indicator of the order of the diagram. 

We are now in a position to develop two tree level quantities - namely the clas- 
sical density and the classical response function. Diagrammatically, we represent 
the classical densities by wavy lines and the classical response functions by thick 
lines. The tree level density (0) is given by the sum of all tree diagrams which end 
with a G^^ propagator, as shown in figure 4(b). This is equivalent to the mean 
field equation, as may be seen by acting on both sides of the graphical equation 
by the inverse Green function 2((?f — V^). Similarly, acting on the much simpler 
tree level diagram for (■0) (figure 4(a)) with the inverse Green function gives its 
classical equation. 

We now define the three response functions for the theory: 

(V'(x, -k^, -s)i^{x\ k±, s))(^) {(l){x, -k±, -s)^{x', k±, s))(^) 

{(pix,-k±,-s)i>{x',k^,s))^'\ (76) 

where the superscript '1' indicates that they are defined in Field Theory I. Their dia- 
grammatic sums are shown in figure 5. The first one: {^{x, —k±, —s)ip{x', k±, s))*-^-* 
is simply the propagator G^^j^, where the superscript '2' indicates that it belongs 
to the second field theory. It is also easy to show that the second response function 
(0(x, —k±, —s)(j){x', k±, s))*^^) is equivalent to the propagator G^^^^. To do this, we 
rearrange the unrescaled equation for G^^^^: 

2(5 + kl- dl)G^^I^{x, x', k^, s) = -4A0e(a;)G^'^^(x, x', k^, s) + S{x - x'). (77) 

Including a delta function integration in the first term on the right hand side, and 



15 



acting on both sides with the inverse Green function, we obtain 

where G^^j is the propagator for the first field theory. Iteration now generates 
the appropriate tree level expansion for (0(a;, — — s)0(a;', s))*^"'^^ and we have 
^If/^ ~ {'P'PY^^ ■ The remaining response function {(f)tp)^^^ is, as would be expected, 
equivalent in the second field theory to the two point vertex sandwiched between 
a (p propagator and a ip propagator. 



4 Density and Reaction Front Calculations 

We first note that we cannot draw diagrams which terminate with a ipi propagator 
in Field Theory II. Consequently, we conclude that (ipi) = 0, and hence that 
{ip) = ipc- This also follows from averaging equation (|37|). We now turn to the 
asymptotic evaluation of (0). Inserting the classical (tree level) solution ( ^If ) into 
the Callan-Symanzik solution ([TTD, and making the leading order replacements 
Qri^'^i and for large qr — ^ g*^, we obtain 

(0) = 1 j|^| + o.3787^|x|(^-^'^)/i2g-(2/%«^^^^"/^l-l^^^^'^^ + . . . (79) 
If we use the explicit value of g*^ from (|66|), then 

= 1 J|^| + o.3787-/^|x|(^-^'^)/i2e-(2/=')(^-^)'^'''"''l-l*''"^''' + . . . (80) 
2 (47re)^'-'-^ 

So the tree level expression consists of the expected linear term, which must be 
present if the boundary conditions are to be satisfied, together with a stretched 
exponential component. 

4.1 One Loop Contributions 

According to our earlier arguments we expect the next order contributions to (0) 
(in Field Theory I) to contain one loop embedded somewhere in the tree diagram. 
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The diagrams corresponding to this prescription are shown in figure 6. However, 
we have shown that we may translate these diagrams into Field Theory 11 by re- 
placing response functions by propagators. This is convenient as we have analytic 
expressions for the Green functions in the second field theory (at least asymptoti- 
cally), and so performing calculations becomes easier. The equivalent diagrams for 
Field Theory II are shown in figure 7. Notice that the density fines present in the 
diagrams for Field Theory I have been absorbed into the vertices for Field Theory 
II, where a factor of 4>l — is present at the source vertex. 

We begin by calculating the loop contained in the third diagram of figure 7 
(whilst not including its leftmost vertex). At this level of approximation, we replace 
the source (the incoming classical density lines at the rightmost vertex in Field 
Theory I) by a delta function at the origin, with a weight equal to the area under 
the classical reaction front; in other words: 



This will be vafid provided (0) decays much more slowly than the classical reaction 
front - an assumption that will be shown to be justified a posteriori. By integrating 
the classical equations, we also have the relation 



which is simply saying that, classically, the number of particles entering the system 
is the same as the number being annihilated at the reaction front. This relationship 
is also true non-classically, if we average over the noise. After we have performed 
the rescaling x' — (Xjy^^x', this becomes 




(81) 




(82) 




(83) 



So the vertex factor at the source becomes: 




(84) 
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Hence the loop is given by the integral 



2\-^'^J^'\\JYI^ I 



-{x'+kl+isy/'^\x'-x'{\ 



4(fc2 +i5)V2 4(A;i-i5)i/2 4(5/ + A;i + z5) V2 



-{x'+ki~isy/^\x'-x'^\ 



(2i'[){2x'i)dx'[dxl 



d'^'^l I ds 



5) 



where the prefactor of '2' counts the number of possible diagram configurations, 
and the s integration is along the real axis. In the integral we have used the form 
of the propagator for the field valid for x'/, Xg, x' ^ 0, and x", x'2 ~ x\ the region 
from which we expect the dominant contribution (as here the propagator falls 
off only as a power law). The x'[ and x'2 integrations are elementary, giving 



2A-i/3p/3(Aj)'^/3 
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g-(fe2^+iS)l/2£'^_(fc2^_js)l/2£/ 



IS 



-1/2 



£'2 



d'^-^k^ds 



(86) 



However, we notice that the leading x' part of the integral is in fact a divergent 
power law. Furthermore, this divergence cannot be cancelled by the renormalisation 
of the theory, as any such cancellation would have to arise from coupling constant 
renormalisation at the tree level (using (|67|)). As the renormalised tree level result 
is still an exponential, cancellation with a power law cannot occur. Consequently, 
we must find another mechanism for the removal of the divergence, and this is 
provided by its cancellation with the divergent loop shown in figure 7b. Turning 
now to the next to leading x' term in the above integral, we have 



_^-l/3J2/3(^^J^d/3^ I ^-(kl+isy/^x' ^-{kl-isy/^i 



1 



+ -r7 



(fci + is)V2 (^2 _ ^5)1/2 



d^-^k^ds 
{27rY 



^7) 



This can be rewritten as 



2 d r 2e-(^l+^^')''"*'-(^i-*^')'^'^' d^-^Lds 



x'^dx'J A{kl + isy/H{kl-isy/^ {2ny 
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i.e. a constant times the derivative of the ip propagator loop integraL Rewriting the 
ip propagators entirely in momentum space, and performing a contour integration 
for s, we end up with 



1 d f e^^^ d'^kdp 
x'^d^ J k2 + (p- k)2 {2ttY+^' 



This integral may be done exactly using some standard results from ||2T]], with the 
result 

j2/3(^j)d/32-i-d^-(rf+i)/2(i _ l^izl^ (90) 

To evaluate the contribution to (0) we now need to include the leftmost vertex and 
propagator: 

- (AJ)^/32-2-d^-(d+l)/2(^ _ I ~,-d-2~-l/2g-xV2|,_,'|^-,^ ^g^^ 

giving the leading order result 

- (AJ)-i2-i-''7r-(^+i)/2(l - d)r {^^^ x"'-^ (92) 

We may now insert this into the Callan-Symanzik solution (|7TD, and use the results 
for the running current /coupling (|7^/ (ffSD, and for the coupling fixed point (|66|). 
This leads to the 1 loop density correction 



^-2d-l 



(93) 



327rVe' 

which justifies the use of the delta function approximation for the source. The one 
remaining diagram in figure 7 consists entirely of (p propagators, and so asymptot- 
ically we expect an exponential dependence which we neglect in comparison with 
the power law. So to this level of accuracy, we have 

W = lj|.| + 0,3787|^|.|<'-«V-e-i(*.)"'^-">l-"- + M^ + ,.. (94) 
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Using the relation (^) it is also straightforward to calculate the form of the reaction 
front profile: 

R = A(0^ - = 0.3787 + i)2|^|(7'^-5)/i2e-§(4-)^/^i-/Vl'-^'/^ 



9/D 

mo 

32M/D 

where only the leading terms generated by the differentiation of each of the com- 



+ (2^±iiP^il^ + .... (95) 



ponent parts of (|9J) have been retained. 

Finally, we consider the cancellation of divergences at 1 loop which, as we 
mentioned earlier, can most easily be seen in the formalism of Field Theory I. We 
expect divergent contributions from the 1 loop diagrams a,b,c and d in figure 6, in 
the limit where the position of the loop's left vertex tends towards that of the right 
vertex. In this limit, where no insertions are possible into the response functions, 
it is appropriate to replace the loop of (p response functions with one of ip response 
functions. Evaluating this loop gives the result l/2g*^, and the diagrams become 
as shown in figure 8. However, if we consider the corrections to the tree level due 
to sub leading terms in goi^gRj (from (|67|) ), we have the same diagrams but with 
opposite signs, which exactly cancel the 1 loop divergences. 

4.1.1 Two Loops 

Whilst we have not calculated in full the contributions to the density from the two 
loop diagrams, a remark concerning their general nature, and of the nature of our 
perturbation expansion, is in order. A sample of these two loop diagrams is shown 
in figure 9. 

The easiest diagrams to evaluate are the first and second of those in figure 9, 
for which it is easy to check that they have the form 

|^|-4d-3 J-3 



e2 



~ e-i/Vn^l^l'+'e^^T'- (96) 
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Hence the perturbative expansion for the power law contributions to (0) would 
appear to have the form: 

{(f)) ~ e^'^lxrVl^f+'e'/^] + e-'^\x\-''[J\xf+h'/'']-' + e-'^\x\-''[J\xf+V/^]~^ + ... 

(97) 

Consequently, we see that the condition for our field theory to be valid is that the 
dimensionless parameter J|x|'^"'"^ be ^ 1. It should also be noted that subleading 
power laws from the loop integrals will be smaller than their leading term by factors 
of (J|a;|^+i)-i/3. 

4.2 d>2 

At the upper critical dimension for the system, in this case d = 2, we expect 
logarithmic corrections to the d > 2 results, owing to the presence of the marginally 
irrelevant parameter A. The Callan-Symanzik solution ([7T|) is still valid, although 
with a different coupling, which we calculate by taking e ^ in (|65|) . This gives 
the running coupling 

gjiU-^) = r^^^— . (98) 

The behaviour of the running current is as previously calculated. Using the asymp- 
totic form g^r^ Atc/ ln(Ka;), we obtain 

_5_ 

f <\ 1 Ti I /ln|a;|\^^, , i 2^. 7u/2n„ uh-i/2|„|3/2 |x|"^ln|x| 

(99) 

where higher order corrections will be only 0[(ln |a;|)~^] smaller, so the asymptotic 
regime will be accordingly hard to reach. Finally, for the reaction front, we have 

^-^■'^^^^^ (ln|x|)Vi2 ^ ' + I6vr2 

(100) 

For dimensions higher than the critical dimension, the expressions from the 
evaluation of the Feynman diagrams are used directly without being inserted into 
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the Callan-Symanzik solution. This gives us the results, valid for d > 2, and in the 
regime {J^X/ D^)\x\'^+^ > 1: 

-(A J)-i2-i-'^7r-("+i)/2(^ _ |^|-d-3 ^ ^ ^ ^ ^^Qi) 

and 

i? = 0.3787p9/i2)^^/i2D|x|=^/^e-§(^^")'''l^l''' (102) 
-DA-iJ-i2-i-'^7r-('^+i)/2r (i_rf)(rf + 3)(rf + 4)|a;|-'^-5 + ___ 



5 Interface Fluctuations 

We now turn to the related problem of the nature of fluctuations in position of the 
reaction front. This is similar to the question of the fluctuations of an interface in 
the dynamical Ising Model, as described by the time dependent Landau-Ginzburg 
(TDLG) equation with noise (for example in model A - see [0). This equation 
may be mapped to a path integral for the field $, with the introduction of response 
fields $, using the Martin-Siggia-Rose formalism: 

where the last term in the action results from averaging over the noise. Solving 
the TDLG equation in the absence of noise gives us the classical profile $c, and on 
physical grounds we expect the full functional form of $ to be — f{y,t)) ^ 
$c(a;) — f{y,t)^'dx). The idea now is to substitute this into the action and to 
expand the response fields in terms of some complete set of eigenfunctions \l/n(a;) : 

^(x, y,t)=Y^ Ajn{y, t)^„(x), (104) 

n 

where the {An} are normalising constants. This set is chosen such that when the 
X dependence is integrated out of the action, it leaves behind an unambiguous 
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equation for f{y,t), obtained by integrating over the new response fields f{y,t) 
in the path integral For the Ising case, f{y,t) can be shown to satisfy a noisy 
diffusion equation, whose solution implies that fluctuations delocalise the interface 
for d < 3. A similar result for reaction fronts would have dramatic consequences. 

Returning to the reaction-diffusion system, we expect the functional forms for 
the fields in our geometry to be ip{x — fi{y, t)) and — f2{y, t)), by analogy with 
the Ising case. Considering first the situation where we neglect noise in the system, 
we expand the above functional forms, giving 

<f> ^ <P{x) - f2iy, t)ct>\x) ^ct>= ^-J\x\ T \jf2 {x = ±L) (105) 
= i^{x) - My, tW{x) ^ V' = -\jx + \jh. (106) 

Hence a = (l/2)J(/i - /a) and b ^ Jx - (l/2)J(/i + /a) at x = L, and a = 
—Jx+ (l/2)J(/i +/2) and b = (l/2)J(/2 — /i) at a; = —L. In the absence of noise 
a and b represent the (positive) particle densities, so we must have fi — f2- 

However, if we include the noise term then this argument is invalid, and we 
proceed, as in the Ising case, by inserting the expanded functional forms for 4> and 
ijj into the action for Field Theory I, giving 

S^l dxd''-'ydt[24>m-f2 + Vi/2] + 2AV'cV'c(/i - /2)} 

+2^K[-/i + Vi/i]} + A(0^ - ^2)(0^ - V.2)], (107) 

where we have made the approximation (p ^ (p^, and then used the classical equa- 
tions to simplify the expression. Here y are the coordinates for directions perpen- 
dicular to the applied currents. In our case it is now appropriate to Fourier expand 
the -0 and fields, i.e. 

= E Uy, t)9n{x) - E Uv. t)0n{x), (108) 

n n 

where 9n — sin(n7rx/L) for n > 0, 9n — 1/2 for n — 0, and 9n — cos{mrx/L) for 
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n < 0. Inserting this into the noise term and performing the x integration, we have 

/ V^2(02 _ ^2^^^ _ ^ ^^^^ Orfimicfl - i^l)dx = J] Cn^^^Cm, (109) 



n,m ^ n,m 



where Xnm is a symmetric matrix which we now diagonahse. Using Cn = DnmCm, 
but such that Co = Co, we rewrite ( |109| ) as 

J2CnKmCm, (HO) 
n,m 

where A is a diagonal, and D a diagonahsing, matrix. Bearing in mind the symme- 
tries of the classical solutions, we can perform the x integration within the action 
to arrive at the path integral: 

/ Vf.VhYl'^CnT^Uexp \- [ d''''ydt\2j2UM-f2 + ^lf2) (HI) 

n L I n>0 

+2XBn{fl — /2)] + X/ ^n^mXnm + 2CoC*o[— /l + V^/l] " A ^ CnA„„C„ 
n,m n j 

with 

An = e^<p'Jx, Bn = Onipc^'Jx, Co = Oo^p'^dx ~ JL. (112) 
Integrating over Co and /i gives the equation 

- /i + Vl/i + r/ = 0, (113) 
where is a (possibly imaginary) noise variable, with a Gaussian distribution: 

P{7]) ~ e-('''^oA|Aoo|)_ ^114) 

If we also diagonalise the noise term involving the ^ fields, then the relevant part 
of the action is transformed to 



2 E WnL[An{9 - Vlg - 77) - 2\B^g] + \Y.LKnL (115) 

n>0 n 



where g = /i — /2, and the equation for fi (|113|) has been added into the action. The 



are coefficients generated by writing C,n in terms of {^m}- Finally, performing 
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the integrations over ^„ and g, we find equations for g which can only be mutually 
consistent for different n if (7 = 0, or in other words, if /i = /2 = /• From (|109D we 
see that Aqo ~ JX^^, and so 



-/ + Vl/ + r/ = 0, 



where 77 is a Gaussian noise variable with probability distribution: 



(116) 



P{r]) ~ e 



-(const.)L'^ Jif' 



(117) 



We now proceed to calculate the mean square fiuctuation (/^) — (/)^. This can be 
done in a straightforward manner, solving the noisy diffusion equation satisfied by 
/ using a Green function method. The results are 



ln(LxA)J~ 



for > 4 
for d = 2) 
for = 2, 



118) 



■d-l 



and A is now the large k 



where the system has physical dimensions L\\ x , 
momentum cut-off. So we expect that interface fluctuations will be unimportant if 

L2(JD-1)1/2 



{f) - if)' 



In(LxA) 
L'iJD 



LiJD 



A('i-3)/2 

L2(JD-1)1/2 

(ln(LxA))i/2 
L|(JD-i)i/2 
7172 



> 1 for d > 4 

> 1 for d = 3 (119) 

> 1 for = 2. 



These results can now be applied to the problem of the late time behaviour of 
an initially homogeneous distribution of A and B particles [23-26], where it has 
been shown that the reactants segregate asymptotically ||23|, |2^. We assume here 
that we can access the quasistatic time dependent regime by simply replacing our 
currents J by their time dependent analogues (this point is discussed further in the 
next section). In |2^ it is demonstrated that these time dependent inward currents 
(towards the domain interfaces) scale as J ~ ^-{'^+2)74^ where the domains have 
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a characteristic length scale which grows in time as t^^^. So on the basis of our 



assumption we can insert the appropriate time dependencies into ( |119| ), from where 
it is easily seen that fluctuations are unimportant for large enough t, in dimensions 
where segregation occurs {d < 4). 

6 Discussion 

The main results of our earlier calculations are expansions for the asymptotic be- 
haviour of the density and reaction front proflles for dimensions above, below, and 
equal to the critical dimension. We now compare our analytic results with the 
available data from recent numerical simulations 0, |^, [1^, |1^ |15|. Note that in 
all of these papers except ||T^, the initial conditions are those of complete parti- 
cle segregation - so the particle currents at later times are time dependent. The 
remaining reference [O contains the results of simulations in the steady state. In 



principle, the calculations of this paper can be redone for the time dependent case, 
but simple one loop considerations for {■ip'^) indicate that the dominant contribu- 
tions to the integrals originate from large times. At these times the reaction front 
is formed quasistatically, and so we expect to be able to relate to the steady state 
case by making the correspondence J ~ t^^/^ [[Tsl , p6| (but see below for occasions 



where this breaks down). Data for d = 2 in the time dependent situation is pre- 
sented in 10] and , although in there is insufficient information to extract the 
asymptotic behaviour of the reaction front. Further simulations for d = 2 and also 
for d = 1,3 are given in |]T^, where evidence for (j^) - their proposed scaling form of 
R- is given. The reaction front profile is seen to exhibit good scaling collapse close 
to its centre for d = 1,2,3 but again no information is available for the asymptotics 
addressed in this paper. 



Turning now to the Id case, the simulations in [14, 15] were performed using 
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an infinite reaction rate constant, i.e. if two particles of different species eitlier 
crossed or occupied tlie same lattice site, they immediately annihilated. With 
initial conditions of complete particle segregation, this resulted in total separation 
of the two species at all later times. Consequently, the reaction front profile was 
determined by the fluctuations in position of a delta function like reaction front. 
Our results are for finite reaction rates, and are dominated by density fluctuations 
which propagate out from the reaction front centre to positions far away, a process 
which cannot occur in the Id simulations mentioned above. We believe this to be 
the reason for the discrepancy between our analytic calculations and the numerical 
results. For example, the Id simulations of Cornell produce evidence for a Gaussian 
reaction front profile, most notably in figure 8 of |T^. In that graph logi? is 
plotted against (x/X^^^)^, where X^"^^ is the width of the reaction product profile 
C = J R{x,t)dt, as measured by its second spatial moment. The resulting straight 
line indicates that the Gaussian profile is maintained well into the asymptotic region 
(i.e. at least as far as {x/X^'^^y ^ 30) - in exactly the region where, in our model, 
we would expect our asymptotic expansion to begin to apply. 

In addition, controversy still exists over the spatial moments of the reaction front 



profile - Araujo et al and Cornell fl^ disagree over the presence of multiscaling. 



In fact, our calculations suggest that multiscaling does indeed occur for high enough 
moments in the time dependent version of our model, starting from completely 
segregated initial conditions. For the steady state situation, the existence of the 
asymptotic power laws found above implies that the moments: 

,(,..f»^)''' (120) 

do not exist for g > + 1. However, in the time dependent case, these moments 
must exist due to the presence of a diffusive cutoff at a;/? ~ t^/^ Therefore, 
for the calculation of the spatial moments x'-^-' (for large enough g), we cannot 
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relate the steady state case to the time dependent case by simply applying the 
scaling substitution J ~ t^^/^. We can make these remarks more quantitative by 
performing the calculation of the spatial moments in the time dependent situation. 
Separate arguments must be applied for d < 2, when our RG arguments imply that 
the steady state profile has a scaling form (|95|) ; and for d > 2, when ( |102|) shows 
that the fluctuation induced power law tails do not scale. For d < 2, we have: 



where a and /3 are defined in the usual way ^ Here F{y) is a function which 
provides a cutoff at y ~ 0(1), but whose inclusion does not affect the calculation 
of moments which are finite even in the absence of a cutoff. For q < u + 1, where 
the qth moment of S is finite (even without a cutoff), we can therefore neglect the 
effects of F. However, for g > + 1, the gth moment is infinite without the cutoff, 
so the integral will now be dominated by the region (x/t^^^) ~ 0(1), where the 
asymptotic result S ~ {x/t°')~^~'^ may be used. These considerations lead to the 

result x^'^^t) ~ t"''', where (neglecting any logarithmic corrections for q = cu + 1): 

{a for g < cij + 1 

1 , (a;+l)(a-l/2) r ^ , i ^^'^'^^ 

i _l_ A ix. — for g > c<j + 1. 

Hence we have a cusp at g = + 1, above which tends towards 1/2 for large 

g. Note that this value of 1/2 is specific to a diffusive cutoff of the form F(x/t^^'^). 

For d = 2 we also expect logarithmic corrections to the above power laws. 

For d > 2, we must carry out a slightly different calculation, as although the 

classical (tree level) reaction front obeys scaling, ( |102| ) reveals that the one loop 

power law correction does not. However, for moments which exist without a cutoff, 

it turns out that the classical terms still give the dominant contribution in the 

scaling limit. For these terms we have, in the steady state case: 

/OO roo 
\x\m,ix, A, J)dx ~ A^/V^/^' / \x\'^Sc[{\JY^^x]dx ~ A-'?/V^-«/^ (123) 
-oo J —oo 
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However, for the non-scaling power law we must consider: 

■/(Aj2)-l/{d+4) ^ ^ 

where we have imposed a lower cutoff in the integral, derived from the expansion 
parameter of the d > 2 asymptotic series (|102|) . Comparing the J dependence of 



the two results above, we see that the first of these will dominate in the scaling limit 
J — > 0. Substituting J ~ t^^/^ and normalising, we end up with Uq = a = 1/6, 
for g < (i + 4. For the higher moments (g > c? + 4) we need to introduce the cutoff 
function F, so the integral will be dominated by the region (x/t^^'^) ~ 0(1), where 
we can use the asymptotic power law from ( |102|) : 

( \f-l/2\-l ,.oo rf-q 

^ / dx4-,n^/t"'') ~ A-H(^-2-'^)/2. (125) 

t~^l^ Jo 

Consequently, we have the result where (neglecting logarithmic correc- 

tions for g = c? -|- 4): 

{I for g < d + 4 

(126) 
i-^ forg>rf + 4. 

In this case we have a discontinuity at g = d -|- 4, a result of the power law term 
being unimportant for q < d + A, but dominant for g > (i + 4. Once again we stress 
that the limiting behaviour ^ 1/2 as g oo is dependent on the diffusive form 
of the cutoff. 

Thus, we predict the existence of multiscaling in the time dependent case in 
qualitative agreement with Araujo et al, even though we are considering a different 
model. In general, power law tails in the steady state reaction front profiles should 
always lead to dynamic multiscaling, whatever their origin. These arguments are 
similar to those of Cornell et al. who find evidence for multiscaling in the 

reaction nA + niB with {n,m) ^ (1, 1). However, in that case the solutions of 
the mean field rate equations already give power laws, even without the addition 
of fluctuation effects. 
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Finally, we conclude that the available simulations are not directly applicable 
to our calculations of asymptotic power laws and multiscaling. However, if the 
asymptotics could be reached in a model with a finite reaction rate, our results 
should be amenable to numerical tests. These might be easiest in Id where the 
power law tail should be most pronounced. 
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2 Vertices for Field Theory II. The couplings associated with each 



of the above diagrams are: (a). 
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Figure 1: Vertices for Field Theory I. 

> '> ^- >- >- > 

(a) (b) (c) (d) (e) (f) (g) 




Figure 2: Vertices for Field Theory II. The couplings associated with each of the 
above diagrams are: (a). -A^^ j4/3/i(^) (b). >}/^J^'^h{x) (c). -2{XJY/^x 
(d). {XJf'^'x (e). -(AJ)2/3x (f). -2(AJ)2/3/(x) (g). 2(AJ)2/3/(x) (h). 
-2A (i). 2A (j). -A (k). A (1). A (m). -A. 



X(k,t) 




Figure 3: The sum of diagrams contributing to the primitively divergent vertex 
function A(A;, t). 
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Figure 4: Tree level diagrams for (a) (■?/'), and (b) (0). 
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Figure 5: Response functions for Field Theory I: (a) {ipip)^^\ (b) and (c) 
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Figure 6: One loop diagrams in Field Theory I. 
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Figure 7: One loop diagrams in Field Theory II. 




Figure 8: Divergences at one loop. The factor underneath each diagram is associ- 
ated with the vertex. 
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